function plot_function_1d(u, a, b)

if nargin == 2
plot(a, u(a));
end

if nargin == 3
x = a:0.01:b;
plot(x, u(x));
end

end